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ABSTRACT 

I 

A finite  element  analysis  was  made  of  crack  growth  in  a 
center-cracked  specimen  subjected  to  monotonically  increasing 
load  until  the  point  of  fast  fracture.  Since  part  of  the 
specimen  experienced  unloading,  the  boundary  value  problem 
which  was  formu3ated  was  based  upon  incremental  theory  of 
plasticity.  lixper imental  load  and  crack-size  records  were 
utilized.  Linear  relations  between  plastic  energy  and  crack 
growth  were  observed.  fracture  toughness  parameters  (>c,  which 
were  evaluated  at  the  onset  of  unstable  crack  propagation,  ob- 
tained from  finite  element  analysis  were  in  good  agreement 
with  those  determined  experimentally. 


1. 


I nt  roduct i on 


Linear  clastic  fracture  mechanics  is  widely  used  by  research 
workers  and  engineers  as  a tool  to  determine  fracture  toughness, 
or  C , for  evaluating  tiic  susceptibility  of  engineering 
materials  to  unstable  fracture.  Since  ductile  materials  have 
higher  resistance  to  crack  growth  than  brittle  materials,  the 
recent  trend  is  to  develop  more  ductile  materials  for  structural 
applications.  However,  ductile  materials  show  a considerable 
amount  of  crack-tip  plasticity  and  significant  amount  of  crack 
growth  prior  to  the  onset  of  unstable  fracture. 

A plastic  zone  correction  [1,2],  the  J-integral  [3,4]  and 
COD  methods  [5]  have  been  proposed  to  treat  fracture  involving 
crack-tip  plasticity  in  the  absence  of  crack  growth.  further, 
the  crack  growth  resistance  curve  method  [6]  has  been  proposed 
for  cases  having  subcritical  crack  growth.  Liebowitz  and 
F.ftis  [7,8]  introduced  the  nonlinear  energy  method  which  is 
applicable  to  scmibrittle  fracture.  Jones  et  al  [9]  developed 
a corresponding  experimental  program  to  determine  the  toughness 
parameter  G^  and  Liebowitz  et  al  [10]  provided  theoretical 
and  experimental  comparisons  between  their  nonlinear  energy 
method  and  other  existing  methods.  Recently  Jones  et  al  [11] 
used  four  empirical  methods  to  determine  (! . in  the  range  of 
crack  growth. 

In  this  work,  a finite  element  analysis  is  developed  for 
crack  growth  problems  based  upon  incremental  theory  of  plasticity. 


In  this  paper,  mention  is  made  of  t lie  work  of  Newman  f 1 2 ] whose 
erack  growth  criterion  is  based  on  crack-tip  strain  with  elenicnt- 
mesh  size  being  a parameter,  and  that  of  do  Koning  [13|  who 
found  crack  tip  opening  angle  as  constant  during  the  process  of 
crack  growth.  In  this  analysis,  the  realistic  experimental 
curve  relating  load  and  crack  size  is  taken  as  an  input.  In 
order  to  demonstrate  the  validity  of  this  work,  the  numerical 
output  of  load  and  displacement  is  compared  to  its  experimental 
counterpart.  Moreover,  a linear  relation  is  observed  between 
plastic  energy  and  crack  size  during  the  process  of  crack 
growth.  The  f in ite- element  values  of  G , based  upon  empirical 
methods  derived  by  Jones  et  al  [11],  are  in  good  agreement 
with  the  experimental  values. 


2.  Stress-Strain  Relations 


In  linear  elasticity,  the  stress-strain  relations  for  a 
homogeneous  isotropic  material  can  he  written  in  either  of  the 
following  forms  [ 1 4 J : 


°ij  = * Ckk  6ij  * 2m  eij  * 


(2.1) 
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(2.2) 


where  X and  p are  the  Lame  constants,  o.^  and  e^j  are  the 
stress  tensor  and  strain  tensor  respectively.  Introducing 
stress  deviator  Sj.  as  follows: 
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i j l j 3 kk  i j 

(2.2)  can  be  rewritten  as: 
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where  H and  v are  Young's  modulus  and  Poisson's  ratio  respec- 
t ively . 

In  nonlinear  elasticity,  adopting  the  model  suggested  by 
Rambcrg  and  Osgood,  we  generalize  the  s t ress - s t ra in  relations 
as  follows: 


lieij  ■ O * V)  si3  • -4-^  °kk  6ij  4 r sij 


e i j , (2.5) 


where  the  effective  stress  oc  is  defined  as: 
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In  simple  tension  test,  namely,  all  stress  components  are  zero 
except  Ojj  =a,  we  obtain  the  following  nonvanishing  components 
of  sjj  and  the  effective  stress: 
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lor  the  theory  of  incremental  plasticity,  a certain  portion  of 
the  material  is  said  to  be  in  an  unloading  situation  if  its  current 
effective  stress  o^  is  less  than  the  maximum  effective  stress 
o*  it  experienced  before,  otherwise  it  is  said  to  be  in  a loading 
situation.  The  relations  between  incremental  stresses  and  in- 
cremental strains  in  the  loading  situation  can  be  derived  from 
eqn.  (2.5)  as  follows: 
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In  the  unloading  situation,  we  assume  the  incremental  stress- 
strain  relations  arc  the  same  as  those  in  linear  elasticity, 
i . c . , 

Wcij  ■ O'vJJsjj  • dokk  6m  . (2.9) 

The  st ress - st ra in  relation  in  the  case  of  simple  tension  is 
illustrated  graphically  in  Fig.  1.  In  the  case  of  generalized 


plane  stress,  cqn.  (2.8)  may  be  written  in  the  fallowing  matri. 
form: 
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Hquivalcntly , cqn.  (2.10)  can  be  expressed  as 
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where  matrix  [d.jj  is  t lie  inverse  of  matrix  (h^jj.  In  the  case 
of  unloading,  we  have 
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3.  Strain  Energy 

In  the  case  of  a simple  tension  test,  if  the  specimen  is 
first  loaded  monoton ica 1 1 y from  o = 0 to  o = o*  and  then  un- 
loaded monotonical Ly  from  a = o*  to  o = 0,  we  obtain  the  follow- 
ing strain  components: 
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Note  that,  in  evaluating  the  first  and  second  integral  of  eqns . 
(3. 1-3. 2),  the  general  incremental  stress- strain  relations  (2.8) 

/v 

and  (2.9)  arc  utilized  respectively.  = e and  e22  = £33  = e 
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so  obtained  in  eqns.  (3. 1-3. 2)  are  called  plastic  strains  which 
are  the  residuals  after  the  applied  load  is  removed.  Comparing 
eqns.  (3. 1-3. 2)  with  eqn.  (2.7) ^ and  eqn.  (2.7)9,  we  may  rewrite 
cqn.  (2.5)  as 


Reij  = R(ci.i  + eiV  * 


(3.3) 


where  the  elastic  strains  and  plastic  strains  are 
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Now  suppose,  in  a general  case,  the  loading  and  unloading 
history  of  a certain  portion  of  the  material  can  be  divided  into 
two  stages.  In  the  first  stage  the  effective  stress  increases 
monoton  ical  ly  to  o*  which  corresponds  to  o’?,  and  in  the  second 

A 

stage  the  effect  stress  decreases  monotonical ly  to  o?  which  corre- 
sponds to  (Jjj.  Then  the  strain  energy  density  <{>  can  be  obtained  as 
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we  notice  the  strain  energy  density  can  be  separated  into  two 
parts : 
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where  the  clastic  strain  energy  density  4>c  is  a function  of  the 

A 

current  stresses  a ^ ami  the  plastic  strain  energy  density  4^ 
is  a function  of  maximum  effective  stress  a*  . The  irrever- 
sible and  dissipative  nature  of  plastic  strain  energy  density 
is  indicated  in  eqn.  [3,6).  In  two  dimensional  problems,  the  elastic 
and  piast ic- strain  energy  per  unit  thickness  over  an  area  A 
can  be  obtained  as: 


(3.7) 
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Center-Cracked  Specimen 


In  this  work,  we  focus  our  attention  on  a rectangular 
plate  of  length  2£,  width  2w,  and  thickness  B,  with  a centered 
line  crack  of  initial  crack  size  2a^  subjected  to  symmetric 
boundary  conditions  (cf.  Fig.  2).  Therefore  only  the  first 
quadrant  of  the  plate  R = [x,y|0  < x < w,  0 < y < £]  needs  to 
be  analyzed.  experimentally  one  may  obtain  a curve  of  applied 
stress  o vs  crack  size.  A realistic  experimental  curve  relating 
o and  a is  shown  in  Fig.  3.  h'e  are  interested  in  the  investiga- 
tion of  the  process  that  crack  size  slowly  increases  from  a = a^ 
to  a = ac  as  the  applied  stress  increases  from  o - 0 to  a = 

Suppose  we  have  two  adjacent  states:  o = and  a = a^  at 

State  l,o=  0^“^  and  a = a^^  at  State  2.  Correspondingly, 
the  boundary  conditions  may  be  specified  as  (cf.  Fig.  4): 

°x  = 0X>,  = 0 on  Sj  = [x  = w,  0 < y < l]  , (4.1) 

o = o1-1^,  a = 0 on  S = ,[y  = £,  0 < x < w]  , (4.2) 

/ Ay  l 

ux  = 0 , oxy  = 0 on  S-  = [x  = 0 , 0 < y < £]  , (4.3) 

oy  = 0 , oxy  = 0 on  S4  = [y  = 0,  0 < x < a^],  (4.4) 

uy  = o , oxy  = 0 on  S5  2 [y  = 0,  a^1-^  x < w] . (4.5) 

In  other  words,  the  boundary  conditions  associated  with  the 
process  from  State  1 to  State  2 can  be  expressed  as  (cf.  Fig.  5): 
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6u  =*  0 , (So  = 0 011  r,  = [y  = 0,  a1  J<x<w],  (4.11) 

7 A 7 ^ 

where  is  the  stress  of  State  1 distributed  along  y = 0 
between  x < 

If  we  have  the  solutions  of  State  1,  then,  for  any  given 

(2)  ( 7) 

adjacent  state  specified  by  o and  a , the  boundary  con- 
ditions for  the  transition  from  State  1 to  State  2 arc  well  de- 
fined in  eqns.  (4.6-4.11).  The  equations  of  equilibrium  and 
the  strain-displacement  relations  in  the  incremental  form  can 
be  written  as: 


(So..  . = 0 , 

ii,i 


6c . . = T(6u • ■ + 6u • ■ ) . 
ij  2 i.J  J , i 


(4.12) 

(4.13) 


liqns.  (4.6-4.13)  with  eqn.  (2.8)  and/or  cqn.  (2.9)  mathemat  ically 
define  the  crack  growth  problem  as  a boundary  value  problem. 

After  the  boundary  value  problem  is  solved,  the  stresses, 
strains,  and  displacements  at  State  2 may  be  obtained  as 


S.  Finite  I. lenient  Procedure 


Let  a finite  element  mesh  with  N noJal  points  and 

P ® 

triangular  elements  be  set  up  in  region  R = [x  ,y  1 0 < w , 0 < y < 2.]  . 

In  this  work,  N = 300,  N = 51S,  and  a series  of  very  fine 
p e 

elements  are  being  distributed  along  the  path  of  slow  crack 
growth.  The  area  of  those  elements  is  in  the  order  of  10  ^wJl. 

If  j is  the  number  of  a certain  nodal  point,  then  u^jj  and  u^j 
arc  the  displacement  of  that  point  in  x and  y direction  respec- 
tively, f...  . and  f^.  are  the  corresponding  external  concentrated 

4 J ’ 1 L J 

force  components  acting  on  that  point.  The  displacement  field 
within  each  triangular  element  is  assumed  to  be  linear  with 
respect  to  the  coordinates.  Therefore  the  strain  field,  and  accord 
ingly  the  stress  field,  within  each  triangular  element  are  constant 
For  each  element  (cf.  Fig.  b)  , let  the  incremental  strain  field 
[fie],  incremental  stress  field  [fio],  and  incremental  nodal  point 
displacements  [6]  be  represented  by 
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[ficl  = [ fi  £ x , 5cy,  fiyxy]  » (5.1) 
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Then  we  have 

[6c]  = [B]  [fi ] , (5.4) 
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where  [ d J is  a 3 x 3 matrix  as  indicated  in  eqns.  (2.11)  and  (2.12) 
which  are  valid  for  loading  and  unloading  respectively,  and 


b. 

i 

0 

b. 

J 

IB] 

II 

sr 

0 

c • 

l 

U 

c . 

l 

bi 

c • 

J 

1 

X . 

1 

^i 

2 A 

= dot 

1 

X • 

J 

yj 

1 

xk 

>'k 

(5.5) 


(5.6) 
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with  the  other  coefficients  obtained  by  a cyclic  permutation 
of  subscripts  in  the  order  of  i,  j,  k.  The  stiffness  matrix 
per  unit  thickness  of  this  triangular  element  linking  incre- 
mental forces  and  incremental  displacements  may  be  obtained  as 
[15,1b]: 

T 

[k]  = [13]  [ d ] [ 13  ] A,  (5.8) 


where  A is  the  area  of  the  element.  Finally  the  governing 
equation  can  be  written  as: 
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where  the  2N()  x 2N^  stiffness  matrix  [kj  is  equal  to  the  sum  of 
local  stiffness  matrices  [k]j,  I = 1, NQ  . 

As  matrix  [ J J of  each  element  depends  on  the  current  stresses 
when  that  element  is  in  a loading  situation,  the  matrix  [KJ  also 
depends  on  the  whole  stress  field.  Thus,  an  iteration  process 
has  to  he  taken  to  solve  the  nonl  incar  matrix  equation  (5.9).  The 
iteration  process  could  be  described  as  follows.  Suppose,  in  a step 
by  step  process,  we  have  the  solutions  of  State  1,  (i.e.,  for  each  ele- 
ment we  have  the  current  stresses  o ^ and  the  maximum  effective 
stress  o*) . Then,  for  each  clement,  we  guess  at  the  stresses  o ^ . of 

State  2.  Then  the  current  effective  stress  o and  the  effective 
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stress  o maybe  obtained  respectively  fr 


om 


o = ( a “ ♦ o - o o ♦ 3o  ) 
e x y xy  xv 


(5.10) 
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A L A fa  A A A im 

(a  +o  - o o + 3o  ) 
v x y x y xy' 


(5.11) 


The  corresponding  matrix  (d)  for  each  clement  may  be  found  from  one 
of  the  three  following  cases: 

* ~ * 

( 1 ) o =o  and  o >o 
e c e e 

For  case  (1),  the  clement  is  in  a loading  situation;  matrix  [d]  is 
assumed  to  be  the  average  matrix  corresponding  to  two  states  of 
stresses , i.e., 


dii  = 


i 


(2)  < o*  and  o < o 

c - e c-c 


(5.12) 
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l:or  case  (.2),  the  element  is  in  an  unloading  situation;  matrix 

1<J|  should  be  the  same  as  indicated  in  cqn.  (2.12). 

A‘  * 

(3)  o <o*  and  o > a 
e e e e 

for  case  (3),  the  element  has  experienced  unloading  and  is  being 
loaded  again,  the  matrix  [d]  is  obtained  as 

d i j = f d^fo)  + U-f)  cLjCO)  , (5.13) 

where 

f = (oc  - o*)/(oe  - oe)  , (S.14) 

and  d • j ( 0 ) denotes  the  matrix  corresponding  to  vanishing 
stress  state. 

After  solving  cqn.  (5.9),  one  may  obtain  the  calculated 

stresses  and  its  corresponding  effective  stress  oc  . The 

iteration  process  will  be  continued  until,  for  each  element, 

the  guessed  stresses  and  calculated  stresses  are  approximately  the 

same  and  the  averaged  percentage  difference  between  c*c  and 

o is  below  certain  allowable  values  of  error.  After  the  itera- 
c 

tion  process  is  completed,  it  is  straightforward  to  calculate 
quantities  of  interest  in  fracture  mechanics. 
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0.  Diss  ipatcd  linergy  and  Crack  C.rowth 


The  irreversible  and  dissipative  nature  oi'  plastic  strain 
energy  density  have  been  indicated  in  Sec.  5.  Alter  the  iteration 
process  for  a certain  incremental  step  is  completed,  we  have  the 
correctly  guessed  stresses  o and  the  corresponding  effective 

A 

stress  a of  a new  state  for  each  finite  element.  If  the  newly 

obtained  effective  stress  is  larger  than  the  maximum  effective 

* 

stress  oe  obtained  in  the  previous  states,  one  should  replace 


a by  a , which  will  be  used  as  the  maximum  effective  stress 
e 7 e 

for  the  next  incremental  step.  Now  we  recall  the  expressions  for 
clastic  and  plastic  strain  energy  density  as  follows: 
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We  notice  that  depends  on  the  current  stresses  and  therefore 

demonstrates  its  elastic  nature.  On  the  other  hand,  <J>^  is  a 

* 


monoton  ical ly  increasing  quantity  because  c is  always  increasing, 
this  corresponds  to  the  irreversibility  of  plasticity  and 
the  dissipative  nature  of  plastic  energy.  After  the  finite 
element  analysis  for  each  incremental  step  is  completed,  it  is 
straightforward  to  obtain  the  total  plastic  energy  P (per  unit 
thickness)  for  the  whole  region  R as  follows: 
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P C4>pA)i  . (6.3) 

i=  1 

for  a stationary  crack  problem,  the  crack  size  is  regarded 
as  a given  quantity  and  hence  the  question  is  to  determine  the 
critical  load  at  the  onset  of  slow  crack  growth.  However,  as 
the  crack  starts  to  grow  in  a stable  situation,  one  has  to 
treat  crack  size  as  a variable.  Strictly  speaking,  one  more 
variable  corresponds  to  one  more  governing  equation  for  the 
system.  Before  one  can  establish  a governing  equation  for 
crack  size,  the  relation  between  crack  size  and  other  fracture 
parameter  has  to  be  proposed  and  tested  [12,13].  Motivated  by 
the  Griffith’s  concept  [17,18]  that  surface  energy  is  propor- 
tional to  its  area,  and  the  energy  for  the  new  fracture  surface 
has  to  be  supplied  from  the  released  strain  energy  of  the  speci- 
men, we  decided  to  take  the  energy  approach.  Since  crack  growth 
is  generally  regaraed  as  a irreversible  process  and  crack  size 
is  an  ever  increasing  quantity,  we  propose  that  crack  growth 
is  a function  of  plastic  energy,  i.e., 

a - a0  = f (P)  . (6.4) 

We  mention  that  Broberg  [19]  proposed  a similar  concept  which 
can  be  expressed  as: 

dUp/da  = dD0/da  , (6.5) 

where  is  the  energy  flow  to  the  process  region  and  is 
the  energy  dissipation  in  the  process  region.  However  Broberg’s 


— 


=3 
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engineer iny  approach  is  based  on  the  so  called  approximate 
path- independence  of  J- integral  even  in  the  cases  of  crack 

growth. 

After  we  analyze  a realistic  crack  growth  problem  with  all 
the  input  data  obtained  experimentally,  we  observe  the  following 
f act : 

a - aQ  = 3(P-P0)  (6,5) 

which  means  the  amount  of  crack  growth  is  linearly  proportional 
to  the  increment  of  plastic  energy.  The  normalized  plastic 
energy  per  unit  thickness  p = PU/o“w2  is  plotted  against  normalized 
crack  size  c s a/w  in  fig.  7.  The  correlation  coefficient 
attains  a high  value  of  0.9997.  Eqn.  (6.5)  may  be  rewritten  as 


p = PQ  + ( a - a q ) / 3 

■ cv  VB)  * ® a 


(6. 6) 
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7.  Fracture  Toughness 


Liebowitz  and  Fftis  [7,8]  developed  a nonlinear  energy 
method  to  evaluate  the  toughness  parameter  G from  a single 
load-displacement  record  for  cases  of  no  crack  growth.  The 
method  may  be  briefly  described  as  follows.  First,  let  the 
experiment  load-displacement  record  be  represented  by 

1:  F n 

V * R * *4)  ’ 


(7.1) 


where  v is  the  load  point  displacement,  F is  the  load,  1/M  is 
the  crack  size  dependent  elastic  compliance,  n and  k are  constants. 
Then  for  a given  critical  load  F , one  may  obtain  a dimensionless 
quantity  C representing  the  nonlinearity  as  follows: 

n- 1 


c = l + —'A  r -4 

^ 1 n+l^M 


The  nonlinear  toughness  parameter  G is  obtained  as: 


G = CG 
c c 


(7.2) 


(7.3) 


where  Gc  is  the  corresponding  linear  elastic  fracture  toughness  and 
may  be  obtained  by  using  Irwin's  G-K  relation,  i.e., 


Gc  = K‘/E 


(7.4) 


And  K.  is  determined  using  the  following  polynomial  expression 


for  center-cracked  specimen  [20] 

? 

Kc  - oc/;a  [1  - 0.1  \ * (J)  ] • 


(7.5) 


Jones  ct  al  [11]  further  developed  four  different  methods  to 
account  for  the  effect  of  subcritical  crack  growth.  We  briefly 
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describe  these  four  methods  and  in  each  method  we  compare  the 
experimental  value  of  G.  with  the  numerical  value  of  G^.  obtained 
by  finite -element  analysis. 

Method  1 

The  experimental  part  of  this  method  can  be  separated  into 
two  steps:  (1)  to  determine  C from  the  experimental  load  dis- 
placement record, (.2)  to  calculate  G^  by  eqns.  (7.4-7. 5)  using 
actual  critical  stress  oc  arid  initial  crack  size.  l-'or  illus- 
trative purpose,  we  take  specimen  #11  (2024-13)  as  an  example 
and  the  data  for  that  specimen  arc  listed  below: 


(7.7) 


| 
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Gc  (Exp)  = 521  lb/in.,  G (FE)  = 473  Ib/in. 

The  difference  is  about  10"<>  ; 2°»  comes  from  the  calculation  of 
Uc  due  to  fact  that  eqn.(7.5)  always  gives  higher  value  for  Kc  , 
and  hence  higher  value  for  (7  from  cqn.  (7.4)  and  81  comes 
from  the  calculation  of  C due  to  the  discrepancy  between 
finite-element  and  experimental  data  in  the  load-displacement  record 

Method  2 

The  experimental  part  of  this  method  involves  the  deter- 
mination of  a hypothetical  load-displacement  curve  assuming 
there  is  no  crack  growth.  Eftis  et  al.  [12]  propose  that 

F = Fc  M(a0)/M(ac)  (7.8) 

being  used  as  the  critical  load  corresponding  to  the  hypo- 
thetical load-displacement  curve.  Lee  and  Liebowitz  [22]  de- 
veloped a computer  program  which  can  be  used  to  analyze  a 
center- cracked  specimen  in  the  nonlinear  range  without  crack 
growth,  and  therefore  it  is  straightforward  to  find  F assuming 
critical  displacement  v to  be  the  same  as  obtained  experiment- 
ally. We  find  that 

Gc  (FE)  = C (FE)  Gc  (FE) 

2 

= 1.1325  x 58 1.77  (14.01/12.8) 

= 518  lb/in.  (7.9) 


No  experimental  value  for  Gc  has  been  reported  by  Jones  ct  al 

HI). 


Method  3 


The  experimental  part  of  this  method  involves  the  deter- 
mination of  M(),  n)f  k(),  and  C accordingly,  from  the  portion  of 
load- d i splacement  curve,  which  has  no  crack  growth,  and  the 
value  of  V from  the  following  equation 


M, 


n, 


c j) 
o 


(7.10) 


The  experimental  value  of  Cl  is  reported  to  be  520  lb/in.  and 

the  f i n i t e- c 1 emeu  t value  of  Cl  is  found  to  be  519.7  Ib/in. 

c 


Method  4 

Utilizing  k^  and  n()  obtained  in  Method  5,  we  recall  eqn. 
(.7.1)  and  write  as  follows 


v 

c 


k0CM£) 

c 


(7.11) 


which  is  used  to  determine  M M(a^)  . The  experimental 
part  of  this  method  is  to  determine  C.  based  upon  M , n^,  kg, 

1;  , and  li  . , based  upon  critical  crack  size  a^.  The  experimental 

value  of  (I.  is  reported  to  be  58]  lb/  in.  However,  we  found 

~ » 

that  G^.  (I'li)  = 491.6  lb/in.,  which  we  believe  is  more  reliable. 

In  summary,  the  finite-element  values  of  G for  four  different 
methods  arc  473,  518,  519.7  491.6  comparing  with  experimental 
values  521,  - , 520,  581  respectively. 
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8.  Conclusions 

In  this  work,  we  focused  our  attention  on  the  finite  ele- 
ment analysis  of  crack  growth  in  a center-cracked  specimen 
subjected  to  monotonical ly  increasing  loads  until  fast  fracture 
occurs.  During  the  process  of  crack  extension,  although  the 
applied  load  is  increased  monotonical ly , part  of  the  specimen 
has  experienced  unloading.  Therefore,  in  this  work  the  boundary 
value  problem  formulated  is  based  upon  incremental  theory  of 
plast icity . 

Since,  in  a crack  growth  problem,  crack  size  is  no  longer 
a given  constant  parameter  and  must  be  treated  as  a variable,  a 
governing  equation  for  crack  size  is  needed.  Before  we  can 
formulate  such  a governing  equation,  an  experimental  curve 
relating  load  and  crack  size  is  taken  as  an  input  and  hopefully 
from  the  output  one  may  find  the  relation  between  crack  size 
and  other  fracture  parameters  derivable  from  the  stress  field. 
Because  of  the  irreversible  and  dissipative  nature  of  plasticity 
and  crack  growth,  we  believe  that  crack  size  should  be  related 
to  plastic  energy.  Indeed,  after  the  completion  of  the  analysis, 
we  do  observe  the  fact  that  the  amount  of  crack  growth  is 
linearly  proportional  to  the  increment  of  plastic  energy.  In 
another  study  of  ours  which  is  in  progress,  the  linear  relation, 
will  be  employed  as  an  input  replacing  the  experimental  load- 
crack  size  curve.  If  the  result  yields  a load-crack  size 
curve  which  is  in  agreement  with  the  corresponding  experi- 


crack 


:s 


mental,  arc  then  one  may  claim  a governing  equation  for 
size. 

Another  aspect  which  is  of  concern  is  the  fracture  tough- 
ness parameter  t’.  . . Liobowitz  an«J  his  coworkers  [S-10J  developed 
a nonlinear  energy  method  and  its  corresponding  experimental 
procedure  for  the  determination  of  (',  at  the  onset  of  stable 
crack  growth.  Recently  Jones  et  al  lllj  presented  four  dif- 
ferent empirical  methods  to  account  for  the  effect  of  subcritical 
crack  growth.  These  methods  enable  us  to  determine  the  fracture 

toughness  parameter  (1  at  the  onset  of  unstable  crack  propa- 

c 

gat  ion  based  upon  the  experimental  load- displacement  record. 

Then  the  f i n i t e- e 1 ement  values  which  were  determined  foi  G were 

c 

compared  to  the  corresponding  experimental  values;  satisfactory 
agreement  was  found. 

The  following  points  should  be  stressed.  Inputs  for  the 
computer  program  include  the  material  constants  a and  n,  in  addi- 
tion to  Young's  modulus  ii  and  Poisson's  ratio  v,  determined  from 
the  simple  tension  test  of  unnotched  specimen.  However,  as  a 
matter  of  fact,  a and  n vary  over  a wide  range  when  one  tries  to 
describe  the  experimental  st ress - s t ra in  curve  by  a three  parameter 
expression  and,  from  our  experience,  the  finite  element  computer 
program  is  quite  sensitive  to  the  input  a and  n.  We  believe  this 
is  the  major  reason  for  the  existing  discrepancy  between  the  experi- 
mental load-displacement  curve  and  its  finite-clement  counter- 
pa  rt  , wh  i ch  has  a 5.21) 


difference  in  displacement  at  the  critical 


load.  Although  there  is  a discrepancy,  we  find  that  the  finite 
element  values  of  based  on  the  empirical  methods  have  less 
scatter  than  the  experimental  counterparts.  (lencrally  speaking 
this  work  provides  support  to  the  nonlinear  energy  method  de- 
veloped by  l.iebowitz  and  lift  is  [7|,  lift  Is  and  l.iebowitz  [ 8 j and 
its  corresponding  experimental  procedures  .Jones  et  al  [9], 
l.iebowitz  et  al  (10]  and  .Jones  et  al  (11],  in  determining  the 
fracture  toughness  values  in  the  range  of  stable  crack  growth. 
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